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1.0  PREFACE 


This  is  the  final  report  pursuant  to  DARPA  Order  No.  7090.  This 
report  is  based  on  a  study  that  was  conducted  at  Raytheon  Company, 
Wayland,  Massachusetts  between  November  1989  and  February  1992.  The 
subject  consisted  of  tomographic  mathematical  ideas  applied  to 
radar  detection.  The  study  was  monitored  by  the  Air  Force  Office  of 
Scientific  Research,  Bolling  AFB,  Washington,  DC,  under  Contract 
No.  F49620-89-C-0116. 

The  objective  of  DARPA  Order  No.  7090  was  to  produce  a 
computational  framework  for  incorporating  Radon  Transform 
techniques  into  radar.  The  particular  goal  is  a  novel  radar  system 
that  can  map  electromagnetic  backscatter  as  a  function  of  range  and 
range  rate  (or  Doppler  shift) . 


2 . 0  INTRODUCTION 


The  idea  of  reconstructing  a  two-dimensional  distribution  from 
corresponding  line-integral  projections,  which  are  available  for 
all  aspect  angles  distributed  over  180  degrees,  was  introduced  in 
1917  by  Johann  Radon.  These  mathematical  results  have  since 
inspired  developments 
in  diverse  fields. 

The  best  known  is  in 
clinical  radiology 
where  it  is  referred  to 
as  Computerized 
Tomography  (or  CT) . 

Tomo  in  Greek  means 
slice  (see  fig.  1) .  The 
value  of  CT  is  an 
ability  to  routinely 
produce  x-ray  images  of 
pathologically 
interesting  cross- 
sections  that  are 
normally  inaccessible 
such  as  illustrated  in 
figure  1,  without  destructive  surgery. 

Radar  astronony,  geophisical  probing  ,  and  non-destructive  testing 
at  industrial  laboratories  are  some  other  fields  where  this  idea 
has  also  been  successfully  exploited  [2-5] .  Moreover,  the 
principles  are  not  unfamiliar  to  radar  engineers  who  have  perceived 
that  synthetic  aperture  techniques  (  SAR  )  as  well  as  ISAR  imaging 
are  based  on  projections  similar  to  Radon's  line-integral 
projections  [6,7]. 

Johann  Radon's  mathematical  method  for  reconstructing  two- 
dimensional  distributions  was  recently  adopted  yet  another  time 
[8,9]  .  This  time  the  application  is  to  reconstruct  distributions  of 
radar  backscatter  intensity  as  a  function  of  range  delay  and 
Doppler  shift  (  or  range  and  range-rate) .  The  basis  for  this 
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innovation  in  radar  imaging  technology  is  that  analogous  line- 
integral  projections  can  be  obtained  for  these  distributions  by 
viewing  targets  with  a  matched  filter  radar  system  incorporating  a 
suite  of  linear  FM  pulses  which  vary  in  their  Chirp  slope.  This 
application  stands  alone  among  the  applications  of  tomographic 
methods  since  the  projection  angles  are  controlled  electronically 
in  contrast  with  mechanical  steering  that  is  required  in  all  the 
other  cases. 


3.0  TOMOGRAPHIC  PROJECTIONS  IN  RADAR 


The  procedure  in  modern  radiology  for  imaging  pathologically 
interesting  cross-sections,  based  on  some  distribution  of  x-ray 
attenuation,  denoted  by  g(r,y),  is  founded  on  the  inversion  of  its 
tomographic  projections.  The  classic  inversion  formula,  which  was 
derived  by  Johann  Radon  is 


flr(r,*) 


d..'  dplx',6)  /dx' 
rcos (♦  -  8)  -  x' 


(1) 


where  the  variables  x,  x',  y,  r,  0,  and  y  are 
identified  in  figure  2  and  the  tomographic 
projections  are  denoted  by  p(x',0).  The 
latter  can  be  represented  mathematically  by 
the  following  line-integrations. 


p(x* ,0)  m  f  f  f(x,y)0(x/  -  (x  cos0  -  y  cos0) )  dxdy  (2) 


where  f(x,y)  =  g(r,y)  and  5()  denotes  a  line-type  8-function. 

The  second  equation  is  also  the  description  for  two-dimensional 
Radon  Transforms.  Thus,  this  describes  the  transformation  from  some 
two-dimensional  distribution  to  a  set  of  the  one-dimensional  line- 
integral  projections.  In  medical  imaging,  these  tomographic 
projections,  such  as  they  are  called,  are  comprised  of  discrete 
samples  of  the  total  x-ray  attenuation  across  a  slice  of 
interesting  pathology.  The  classic  configuration  consists  of 
penetrating  x-rays  that  are  parallel  and  equally  spaced  so  that 
there  will  be  appropriate  linear  sampling  of  the  cross-section. 
Regarding  individual  projections,  the  paths  traversed  by  the  x-rays 
are  defined  by  x'=  x  cos0  -  y  sin0,  where  0  denotes  the  angle  of 
penetration;  this  is  held  constant  while  the  individual  projection 
is  being  constructed.  Repeating  such  constructions,  for  a  sequence 
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of  mechanical  realignments  providing  appropriate  sampling  over 
180°,  yields  the  projections  that  are  transformed  into  an  image  of 
the  original  x-ray  attenuation  distribution. 

For  radar  imaging,  tomographic  projections  can  be  obtained  by 
transmitting  a  suite  of  linear  FM  pulses.  To  explain  this  idea, 
the  backscatter  generated  by  one  of  these  pulses  and  intercepted  by 
a  matched  filter,  as  illustrated  in  figure  3,  is  depicted  as  a 
composite  echo  corresponding  to  a  dispersal  of  point  scatterers 
whose  ranges  and  range  rates  are  distributed  in  some  random  way. 


FIGURE  3 

A  RADAR  SCENARIO 
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In  this  particular  model,  the  scattering  scenario  can  be 
represented  by  a  statistical  reduction  of  the  composite  echo 
amplitudes  (C)  distributed  over  range  and  and  range  rate.  Thus, 
expressed  in  terms  of  range  delay  (X)  and  (for  the  narrowband  case) 
in  terms  of  Doppler  shift  ($) ,  an  appropriate  representation  for 
the  scattering  scenario  is  the  scattering  function 


s( x,+)  =  E{  | C(t, t)  |2  }  (3) 


where  E{*  •  •  }  is  utilized  to  denote  the  statistical  expected  value. 
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The  scattering  function  will  also  change  with  time,  but  the 
variations  will  normally  be  relatively  narrowband  with  respect  to 
the  dwell  that  is  required  to  obtain  sufficient  projections.  For 
the  case  of  uniform  individual  scattering  cross-sections  ■  the 
scattering  function  has  all  the  properties  of  a  joint  probabilty 
density* . 

The  scattering  function  is  the  radar  analog  with  respect  to  the 
distribution  of  x-ray  attenuation  in  medical  imaging.  The  following 
provides  the  basis  for  obtaining  tomographic  projections  with  a 
suite  of  linear  FM  pulses. 

The  composite  backscatter  described  above  can  be  expressed 
mathematically  as  follows 

m  m 

C<  t)  -  ff  C(T, 4>)Y(t+T,*)  tfcdt  (4) 


where  'F(t,$)  denotes  the  complex  notation  for  transmission  signals 
displaced  with  respect  to  the  carrier  f0  by  $  Doppler  shift.  The 
classic  criteria  for  the  narrowband  approximation  of  the  Doppler 
effect  is  that  time -bandwidth  product  must  be  smaller  than  half  the 
velocity  ratio  c/v  where  c  and  v  denote  light  and  range-rate 
velocities,  respectively.  This  condition  is  almost  always  met  by 
radar  systems  in  operation  today. 

A  matched  filter  that  conforms  with  respect  to  the  description  of 
the  transmitted  waveform  will  yield  the  following  response 
corresponding  to  (4) . 

m 

Y(t>  -  /  C (v) Y* ( t+v)  dv 

-m 

m  m  m 

«  J* J j  C(t,$)  Y(t+v,$)  Y*(t+v)  dvdxdfy  (5) 


By  substituting  v+x=e,  (5)  can  be  rewritten  as  follows. 

m  m  m 

y(t)  =  fjc(i.t)  [f  Y(e,4»T*(t -T+t)<fc]  did*  (6) 


If  n(t)  is  the  complex  pre-envelope  modulation  superimposed  on 
transmissions  of  a  carrier  f0,  then 

¥(«,♦)  =  i»(e)  exp  [j2%  (f0+$)e] 

Y(t-x-e)  *  ji(t-x+e)  exp[i2x(f0(t-x+€)] 
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An  equivalent  mathematical  representation  for  (6)  can  now  be 
derived  by  introducing  Woodward's  two-dimensional  correlation 
function. 


*  j*  n(t)|»*(tn)  exp(72x$t>  dt  (7) 


Thus,  it  is  possible  to  finally  express  the  matched  filter  output 
as  follows. 


Y(t)  *  JJ  C(t,$)  ex^-j2itf0(x-t))  did}  (8) 


Based  on  the  stochastic  model  that  was  adopted  above  for  the 
backscatter  amplitudes,  this  matched  filter  output  represents  a 
stochastic  process.  This  can  be  reduced  by  invoking  the  following 
ensemble  averaging. 


-  //  |x<T-e,*>|*  *<#  (9) 


where  |x(?,$) I2  denotes  the  radar  ambiguity  function. 

The  similarity  of  equation  (9)  with  medical -imaging  projections  can 
be  revealed  by  examining  the  ambiguity  function  for  linear  FM 
(LFM) .  The  characteristic  ambiguity  function,  with  respect  to 
linear  FM  pulses,  has  the  appearance  of  a  ridge-like  surface  with 
an  alignment  which  coincides  with  vt-$=0  where  v,  the  slope  of  this 
line,  is  equal  to  the  LFM  (  or,  Chirp  )  rate.  (Note  that  the  symbol 
v  was  previously  used  in  this  document  with  a  different  identity.) 
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The  following  ambiguity  function  properties  provide  another 
perspective  of  the  relationship  between  the  direction  of  this  ridge 
and  the  Chirp  rate. 

PROPERTY  I 

IF:  |ftU>  "  x(*»*> 

THEN:  |i(t)  exp(jfat2}  -  X  (*»♦"—) 

ic 


IF: 

TH^N: 


PROPERTY  II 
U(f)  - 

U(f)  expljhf2)  - 


&  H(t)  ) 
x(t+**,4>> 

1C 


Figure  4  provides  a  specific  example  of  the  ambiguity  function 
associated  with  Chirp  pulse  compression,  computed  for  a  time 
duration  (T) -swept  bandwidth  (Af)  product  of  15.  Only  the  upper 
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right  quadrant  of  (11)  is  displayed.  Since  the  base  coordinates  X 
and  <|>  are  normalized  in  this  figure,  with  respect  T  and  Af, 
respectively,  the  ridge  line  may  appear  to  have  unity  slope  but 
actually  the  slope  can  be  anything  depending  on  how  the  product  TAf 
is  factored  into  a  ratio. 

For  pulses  with  sufficiently  large  bandwidth,  the  ridge  will  be 
sufficiently  sharp  that  lx(X,$) I2  can  be  approximated  by  the 
line-type  8-function,  so  that  (9)  reduces  to  the  output 


fi t,p)  «  J f  E{ |C(t ,$)  |a} 8  ( t—  (x+-^ ) )  dx d* 


(12) 


where  P=arctan{v} .  Some  interesting  projection  options  are 
presented  in  Section  6  of  this  report . 


The  illustration  in  figure  5.  will  be 
helpful  regarding  an  interpretation  of 
(12)  .  The  ambiguity  function  is  depicted 
in  this  figure  as  a  v-tilted  window  that 
travels  laterally  across  the  x<|>-plane 
parallel  to  the  x  axis.  The  total 
scattering  power  overlayed  by  this  window 
is  effectively  accumulated  at  every 
lateral  position  as  a  function  of  time. 
Thus,  this  output  may  be  interpreted  as  a 
tomographic  projection  of  the  scattering 
function,  therefore  it  can  be  combined 
with  other  such  projections  according  to 
the  inversion  formula  given  by  equation  1. 
However,  to  account  for  the  relationship 
t=x'/sin  P,  equation  1  must  be  amended  as 


Contour  of  radtr  echo 
intensity 


FIGURE  5.  RADAR - ANALOGY  PROJ1CTIORS 


follows . 


M 

2  • 


gU.*)  -  -ji;/  <K>/ 

*7 


dt 


8f(  t,p)/8t 

rsin(^r-P)  -t  sinP 


(13) 


where  g  (r, y)  =s  ( x ,  <|»)  . 

The  projections  that  are  combined  according  to  (13)  can  be  obtained 
by  transmitting  a  suite  of  linear  FM  pulses  which  vary  in  their 
Chirp  rate;  the  suite  size  is  one  of  the  issues  discussed  in  the 
next  section.  Since  different  perspectives  of  the  scattering 
function  are  provided  by  changing  the  Chirp  rate,  this  application 
of  tomography  differs  from  other  applications  in  which  mechnical 
rotation  techniques  are  necessary. 

Commenting  on  the  stochstic  backscatter  model  that  was  chosen  for 
the  derivation  above,  this  model  was  adopted  because  it  provided  an 
anology  in  radar  for  imaging  with  tomographic  methods.  All  extended 
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targets,  which  consist  of  dense  scattering  points  that  are 
physically  disturbed  by  some  random  process,  are  represented  by 
this  model.  Some  examples  are  (1),  the  ionized  wake  attached  to  a 
rentry  vehicle,  (2),  ionospheric  frequency  versus  altitude 
profiles,  (3)  ,  microburst  turbulence  over  airport  runways,  (4), 
rain  clouds,  (5) ,  radar  clutter  from  foliage,  and  (6) ,  rough¬ 
surfaced  rotating  objects.  The  stochastic  model,  therefore,  is  not 
regarded  as  very  restrictive  in  finding  specific  radar 
applications.  However,  another  model  could  have  been  adopted  which 
possisibly  is  best  characterized  by  the  phrase  "what  you  see  is 
what  you  get".  The  proverbial  point  target,  for  example,  is 
represented  by  this  model.  The  model  can  also  be  applied  to 
represent  several  isolated  scattering  objects.  In  fact,  it  may  be 
applied  for  any  target  scenario  which  is  to  be  imaged  by  snapshot 
measurements;  even  dense  backscattering  falls  in  this  class.  In 
this  case,  the  projections  of  C(t,$)  or  I  C(T,$)I  2  available 
directly  from  the  output  of  the  matched  filter  do  not  represent  a 
stochastic  process  and  so  are  not  reduceable  by  ensemble  averaging. 
A  question  arises:  Is  it  necessary  to  amend  the  inverse  Radon 
transform  for  this  case?  One  way  to  find  out  is  to  engage  in 
controlled  experiments  with  live  C(T,$)  distributions. 
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4.0  INVESTIGATIONS  AND  RESULTS 


To  prove  the  concept  behind  this  novel  application  in  radar  imaging 
of  tomographic  mathematical  ideas,  the  investigation  plan  (figure 
6)  contains  several  avenues  that  can  be  taken  to  produce  and 
process  experimental  tomographic  projections.  This  includes  various 
computational  paths  for  producing  pre-formulated  tomographic 
projections,  two-dimensional  radar  distribution  modeling  as  well  as 
echo  simulations  and  real  backscatter  demonstrations  with  different 
degrees  of  experimental  control. 


The  discussion  in  this  section  concerns  the  conduct  of  this 
investigation  during  the  first  phase  which  consisted  exclusively  of 
computational  procedures  involving  pre- formulated  and  computed 
projections  representing  post  match  filter  responses  obtained  with 
different  Chirp  rates.  Several  issues  are  addressed  and  various 
results  are  documented.  In  particular,  the  following  things  are 
discussed  in  this  section. 

(1)  The  point  spreading  characteristic, 

(2)  Parametric  sensitivities, 

(3)  Adjacent  target  detection, 

(4)  Incomplete  projection  data, 

(5)  The  Radon  transform  approximation. 

During  the  second  phase  in  this  investigation,  which  is  already  in 
stages  of  implementation,  the  plan  is  to  simulate  and  process  pre- 
match  filter  signals  contaminated  by  noise.  The  implementation  for 
this  phase  will  incorporate  provisions,  also,  for  an  investigation 
concerning 

(1)  the  paradoxical  moving  object, 

(2)  coherent  versus  noncoherent  tomographic  imaging. 

The  third  and  the  last  phase  planned  for  this  proof-of-concept 
program  is  test  and  demonstration.  The  conduct  of  this  phase  will 
provide  controlled  experimental  results  for  real  objects.  One  of 
the  system  configurations  described  in  Section  6  can  be  employed  as 
the  basis  for  supporting  hardware.  This  phase  would  have  to  be 
accomplished  as  a  follow-on  of  the  present  work.  Several 
engineering  tasks  are  required,  including 

(1)  system/parameter  trade-offs, 

(2)  test  configuration  trade-offs, 

(3)  hardware- software  developments, 

(4)  test  configuration  integration. 

Initially,  the  demonstrations  can  be  conducted  with  point-like 
object  configurations  consisting  of  one  or  more  targets  of 
opportunity.  The  purpose  for  this  step  is  to  confirm  the  basic 
imaging  capability,  offered  by  tomographic  methods.  The  next  step 
can  be  driven  by  a  more  ambitious  goal  consisting  of  quantifying 
the  ability  to  produce  accurate  range  and  Doppler  estimates.  This 
will  require  greater  control  of  the  test  configuration.  Possible 
approaches  are  (1),  calibration  by  another  radar,  (2),  command- 
control,  and  (3),  a  controllable  platform,  such  as  a  rotating  disc- 
type  surface,  impregnated  with  scatterers.  The  ultimate  step  will 
be  to  demonstrate  high  resolution  imaging  and  to  compare  the 
results  with  conventional  pulse  Doppler  images.  An  expansion  of  the 
follow-on  plan  in  terms  of  greater  detail  is  currently  being 
developed.  Some  preliminary  details  are  provided  in  Appendix  A. 


-10- 


4.1  THE  POINT  SPREADING  CHARACTERISTIC 


The  point-spread  characteristic  and  the  radar  ambiguity  function 
are  interchangeable  in  the  context  of  the  following  discussion. 

The  radar  ambiguity  function  was 
introduced  shortly  after  WW  II  by 
P.M.  Woodward  in  his  famous  monograph 
published  in  1953  [10]  .  Almost  from  the 
start  systems  engineer  were  dedicated 
to  minimizing  detection  ambiguities, 
maximizing  joint  parameter  estimation 
accuracy,  and  maximizing  resolution  and 
quality  imaging  fllj  .  An  endless  search 
began  for  optimum  match-filter 
waveforms  to  employ  in  given 
situations,  and  radar  scientists  were 
forever  intrigued  by  the  thumbtack-like 
ambiguity  characteristics  (  figure  7) . 

Due  to  a  fundemental  constraint  that 
restricts  the  ambiguity  function  volume 
to  equal  normalized  waveform  energy,  the  thumbtack  ambiguity 
function  has  frequently  been  an  ultimate  objective  in  this  search. 

Although  the  waveform  that  can  be  characterized  by  its  thumbtack 
ambiguity  function  has  never  been  found,  certain  large  time- 
bandwidth  waveforms  have  been  discovered  that  qualify 
approximately.  The  shift-register  sequences,  for  example,  supply  a 
code  group  that  has  been  a  popular  source  for  discrete  phase 
modulation,  yielding  an  rms  sidelobe  level  that,  averaged  over  all 
codes  in  the  group,  equals  approximately  one  over  the  square  root 
of  time -bandwidth  product.  Figure  8  displays  the  ambiguity 
functions  that  were  derived  for  two  of  these  waveforms. 
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FIGURE  8.  AMBIGUITY  FUNCTION:  SHIFT  REGISTER  SEQUENCES 
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By  applying  tomographic  methods,  a  composite  ambiguity  function  has 
been  derived  that  possess  the  closest  known  resemblance  so  far  of 
the  elusive  thumbtack.  This  is  presented  in  figure  9  which  was 
derived  by  assuming  a  point  distribution  at  the  origin  and 
performing  the  reconstruction  on  a 
128x128  imaging  palet  with  the 
projections  for  a  suite  of  36  linear  FM 
pulses.  The  projections  consisted  of 
the  discete  samples  of 


P<u) 


2 [ coa  H  ~  l  +  3^nu] 


u 


u 


u  =  nVfx 


(14) 


where  Vf  denotes  bandwidth  (  see 
figure  10)  .  This  function  defines  the 
inpulse  response  for  an  even  symmetric 
ramp  filter  whose  bandwidth  is  Vf.  The 
samples  were  evaluate  at  x=k/Vf, 
k=l,2,3,etc,  so  they  represent  Nyquist 
sampling.  Equivalent  sampling  intervals 
are  employed  in  medical  imaging  during 
conventional  tomographic  processing. 


FIGURE  9.  AMBIGUITY  FUNCTION: 

CHIRP  DIVERSITY  RADAR 
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It  is  possible  to  increase  the  resolution  of  the  point  spread 
function  by  sampling  (14)  at  higher  than  the  Nyquist  rate.  The 
effect  on  the  point  spread  function  is  illustrated  in  figure  11, 
and  several  specific  cases  showing  the  increased  sampling  density 
are  also  presented  in  figure  12.  Positioned  above  each  point  spread 
function  in  figure  12  is  its  corresponding  bandwidth-limited  ( 
Chirp  Diversity  Radar  )  projection;  linear  interpolation  was 
employed  to  interconnect  the  discrete  samples.  The  telescopic-like 
changes  that  are  conspicuous,  as  the  resolution  is  increased,  can 
be  explained  by  noting  that  equal  image  dimensions  were  employed 
for  everyone  of  these  cases. 
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The  analog  form  of  the  point  spread  function  is  revealed  when  the 
resolution  becomes  infinite.  This  is  equivalent  to  Airy's 
celebrated  mathematical  description  for  Fraunhofer  diffraction 
associated  with  circular  apertures. 

fir)  ■  2  i  z  ■  2*Vft  (15) 


where  denotes  the  first  order  bessel  function. 

The  Airy  pattern  is  plotted  in  figure  13,  revealing  potentially 
greater  point  target  resolution  compared  to  the  resolution  that  is 
achieved  with  a  conventional  Chirp  pulse  compression  system.  The 
improvement  factor  is  approximately  2k  based  on  the  mainlobe  widths 
measured  at  the  base.  Another  discovery  is  that  this  composite 
ambiguity  function  is  not  constained  by  uncertainty  principles. 
Hence,  *  the  sands  of  the  beach  ■  is  not  a  suitable  analogy  to 
describe  the  ambiguity  funtion.  Most  important,  however,  the 
composite  ambiguity  function  is  not  affected  by  the  reciprocal 
relationship  between  time  and  frequency  since  the  mainlobe  width  on 
this  is  exclusivly  a  function  of  the  bandwidth.  It  is  conclude, 
therefore,  that  by  employing  tomography  super  resolution  of  radar 
backscatter  is  finally  within  reach. 
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4.2  THE  ALIASING  PHENOMENA 


The  aliasing  phenomena  in  computed  tomography  has  previously  been 
investigated  and  discussed  by  several  authors.  A  list  of  some  of 
these  works  can  be  found  in  the  bibliography  provided  by  S.R.  Dean 
[13]  .  It  was  decide  to  revisit  this  subject  since  it  reveals 
certain  considerations  affecting  the  selection  of  system 
parameters.  It  is  believed  that  the  investigations  discussed  below 
will  provide  some  new  perspectives  on  this  subject. 

Functionally,  the  AIRY  pattern  and  the  (  tomography  )  ambiguity 
function  depart  at  some  point,  depending  on  the  number  of 
projections  that  are  employed  during  tomographic  computa tat ions. 
Figure  14  depicts  three  results  for  which  a  different  number  of 
tomographic  projections  was  employed.  These  tomographic 
constructions  consist  of  64x64  points  and  the  number  of  projections 
corresponding  to  cases  a,b,and  c  consist  of  9, 18, and  36, 
respectively.  (  The  constructions  that  are  displayed  in  figure  14 
were  computed  based  on  the  same  number  of  Nyquist  samples  on  a 
projection  as  the  rectangular  dimension  of  the  construction) . 

Two  distinct  regions  are  visible  in  each  of  the  three  constructions 
that  are  displayed.  The  demarcation  between  these  regions  consists 
of  the  onset  of  ray-like  contour  lines,  an  aliasing  phenomena  that 
has  been  characterized  as  “streaking"  in  the  literature  [14] .  The 


radius  of  the  demarcation  depends  upon  how  many  tomographic 
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projections  are  employed,  and  the  functional  variations  out  to  this 
radius  can  probably  be  shown  to  conform  with  the  AIRY  pattern.  The 
streaks  have  a  relative  amplitude,  compared  to  the  peak,  that 
approximately  equals  the  reciprocal  of  the  projections.  The  step  in 
amplitude  at  the  demarcation  arises  since  the  decreasing  sidelobes 
of  the  AIRY  pattern  are  usually  much  smaller  at  the  demarcation. 

One  puzzling  fact  is  that  the  radius  of  the  demarcation  in  figure 
14c  approximately  equals  one  half  the  rectangular  dimension  of  this 
image.  This  appears  to  contradict  a  prediction  that  was  previously 
reported  in  the  literature.  It  states  that  artifacts  due  to 
aliasing  can  be  avoided  if  the  number  of  projections  is  equal  %/ 2 
times  the  rectangular  dimension  on  which  an  image  is  constructed. 
No  explanation  can  be  given  at  this  time  for  this  discrepancy. 

Some  additional  results  that  were  derived  regarding  this  subject 
are  presented  in  figures  15  and  16.  Figure  15  shows  several 
tomographic  constructions  that  were  derived  with  32,  64,  and  128 
Nyquist  samples  but  18  projections  in  each  case.  Figure  16  depicts, 
for  the  case  of  high  density  sampling  which  was  lOx  the  Nyquist 
rate,  how  tomographic  constructions  are  effected  when  the  number  of 
projections  is  varied.  The  figure  contains  the  tomographic 
constructions  that  were  derived  with  1,  2,  4,  9,  18,  36,  and  72 
projections  employing  128  samples  in  each  case. 


FIGURE  15 

CHIRP  DIVERSITY  RADAR 
POINT  SPREAD  FUNCTIONS 
FOR: 

(a)  32  SAMPLES 

(b)  64  SAMPLES 

(c)  128  SAMPLES 
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The  illustrations  below  provide  additional  examples  of  streaking. 
The  ellipse  is  constructed  on  128x128  points  with  18  projections. 
Both  circle  constructions  also  consist  of  128x128  points.  However 
the  first  circle  was  constructed  with  18  projections  and  the  second 
circle  with  36  projections. 


ADJACENT  TARGET  DETECTIONS 


The  tomographic  images  for  adjacent  point  targets,  denoted  by  A  and 
B,  are  shown  in  figure  17.  These  target  were  positioned 
sufficiently  apart  so  as  to  put  their  images  in  a  neighbor's 
aliasing  artifacts.  The  amplitudes  of  target  B,  corresponding  to 
a,b,  and  c,  respectively,  consist  of  0,  -13DB  and  -25DB.  These 
results  were  derived  by  employing  tomographic  projections  th*fc  w#re 
»pacod  ioq  apart  and  were  comprised  of  64  line-integrations.  The 
relative  amplitude  of  the  streaking  due  to  aliasing  is  slightly 
below  -25DB.  With  no  noise  added  to  interfere  with  the  detections, 
target  B  clearly  is  perceivable  at  -25  DB. 

For  figure  18,  target  B  was  moved  in  closer  to  target  A  so  as  to 
put  both  images  in  regions  unaffected  by  aliasing.  Two  situations 
are  presented  in  figure  18,  each  representing  a  different  distances 
between  the  targets.  In  each  case,  target  B  has  an  amplitude  that 
is  equal  to  -25DB  and  is  competing  with  AIRY  sidelobes. 
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FIGURE  17 

ADJACENT  TARGETS 
WITH  DIFFERENT 
AMPLITUDES: 

(a)  0  DB 

(b)  -10  DB 

(c)  -25  DB 


(C) 
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Figure  19  displays  a  point  target  image  with  certain  projections 
omitted.  In  particular,  it  presents  two  cases  where  projections 
were  omitted  from  a  20°  sector  around  0°-  A  condition  such  as  this 
can  actually  arise  in  practice,  since  it  may  not  always  be 
practical  to  generate  Chirp  pulses  that  will  yield  valid  radar- 
analogous  tomographic  projections  within  narrow  sectors  around  the 
range  delay  x  axis.  Similarly,  but  for  a  different  practical 
reason,  narrow  sectors  around  the  Doppler  axis  $  may  also  present 
a  problem.  In  this  case,  it  may  not  be  practical  to  generate  the 
high  speed  linear  FM  that  is  required  to  provide  projections  near 
90°.  However,  a  short  time  duration,  CW  pulse  is  an  alternative,  in 
lieu  of  an  infinite  linear  FM  rate,  that  can  be  employed  to  obtain 
the  projection  at  90°. 

The  two  cases  presented  in  figure  19  represent  (a)  projections 
separated  by  5°  with  Nyquist  sampling  that  yielded  128  samples,  (b) 
projections  separated  by  10°  with  Nyquist  sampling  that  yielde  64 
samples.  Figure  20  illustrates  what  happens  to  the  image  of  two 
point  targets,  with  one  target  positioned  so  that  its  image  would 
appear  in  the  wedge-like  response  void  which  is  caused  by  missing 
projection  during  tomographic  processing. 
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FIGURE  19 


CHIRP  DIVERSITY  RADAR 
POINT  SPREAD  FUNCTION 
WHEN  PROJECTIONS  ARE 
OMITTED. 


FIGURE  20 

CHIRP  DIVERSITY  RADAR  IMAGE  FOR 
FOR  TWO  POINT  TARGETS,  WITH  ONE 
OF  TARGET  POSITIONED  SO  ITS 
IMAGE  WOULD  APPEAR  IN  THE 
WEDGE-LIKE  RESPONSE  VOID  WHICH 
IS  CAUSED  BY  HISSING  PROJECTIONS 
«  I  DURING  TOMOGRAPHIC  PROCESSING. 


4 . 4  FUZZY  PROJECTIONS 


Exact  resolution  of  radar  backscatter  by  tomographic  methods  is 
achievable  only  with  the  idealized  8- function  model  that  was 
employed  for  the  derivation  of  the  medicine-radar  analogy.  Hence, 
the  tomographic  inversion  formula  is  only  an  approximate  procedure 
for  reconstructing  the  two-dimensional  radar  distributions,  when  it 
combines  the  projections  that  are  actually  generated  by  a  Chirp- 
diversity  radar.  These  projections  can  be  characterized  as  fuzzy 
radar-analogous  projections. 

An  accetable  reconstruction  can  still  be  achieved,  however.  Barring 
practical  constraints,  this  can  be  accomplished  by  employing  a 
swept  frequency  bandwidth  that  is  sufficient  with  respect  to  the 
desired  resolution.  An  alternative,  but  less  direct  approach  is  by 
employing  some  type  of  restoration  that  will  transform  the  fuzzy 
projections  to  more  acurate  tomographic  projections. 

Figures  21  and  22  provide  some  results  concerning  this  issue. 
Figure  21  illustrates  the  point  spread  characteristic  which  was 
constructed  with  fuzzy  projections.  For  this  construction,  the 
projection  corresponding  to  a  point  distribution  was  convolved  with 
(sin{x}/x)2.  Conventional  tomographic  processing  was  then  applied 
The  figure  displays  the  point  spread  characteristic  before  and 
after  affecting  this  transformation.  Widening  of  the  mainlobe  is 
apparent . 


Figure  22  provides  the  tomographic  reconstructions  which  were 
obtained,  with  exact  and  fuzzy  projections,  for  a  hemispheric 
distribution.  Sinct  this  distribution  is  symmetrical,  its 
projections  are  isotropic  and  ideally  described  by  an  ABEL 
transform 


-|*2  ( R 2  -  x2) 


where  R  denotes  its  radius.  For  this  experiment,  the  conventional 
filtering  was  modified.  In  lieu  of  the  conventional  filter,  a 
filter  was  employed  whose  impulse  response  was  equal  to  the 
filtered  projection  from  the  preceding  case.  Effectively,  this 
modification  of  the  filter  had  the  same  effect  as  computing  the 
projections  with  fuzzy-line  integrations.  The  effect  is  barely 
apparent.  This,  probably,  is  due  to  the  nature  of  the  distribution. 
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4.5  A  PRIMATIVE  TARGET  MODEL 


The  rotating  object  is  a  classic  radar  target  that  has  interested 
scientists  as  well  as  research  and  development  engineers  for  many 
years.  Since  target  features  in  crossrange  coincide  uniquely  with 
Doppler  shifts,  the  resolution  of  these  features  is  accomplished 
with  conventional  coherent  Doppler  processing.  The  technique  is 
called  ISAR,  short  for  inverse  synthetic  aperture  radar.  In 
conventional  ISAR  systems,  the  range-Doppler  images  are  derived  by 
illuminating  the  rotating  target  with  a  coherent  sequence  of 
wideband  pulses,  then  processing  the  the  pulsed  echoes  coherently; 
the  generic  implementation  consists  of  a  range-gated  filter  bank  ( 
see  figure  24  )  . 

A  primative  target  model  is  described  in  this  section,  that  has 
been  employed  to  explain  to  others  how  this  imaging  task  can  also 
be  accomplished  by  using  tomographic  methods.  Referring  to  figure 
23a,  the  model  consists  of  a  rotating  ring.  This  is  in  the  same 
plane  as  the  radar  and  visible  to  it  at  every  point  around  the 
cicumference  of  the  ring. 

The  distribution  over  the  -plane  that  this  primative  model 
generates  consists  of  an  elliptic  line-like  distribution.  The 
description  of  this  ellipse  in  terms  range  and  range  rate  is 

(R  ~  r)a  +  <|)a  -  ra 

6 


where  R  and  R  denote  range  and  range  rate,  respectively,  and  r 
denotes  the  radius  of  the  ring,  and  0  is  the  rotation  rate.  This 
distribution  is  drawn  in  figure  23b. 

For  this  distribution,  the  tomographic  projections  can  be  described 
by  a  mathematical  formula  that  reflects  the  dependance  on  the 
particular  viewing  angle  which  is  expressed  as  follows 

t|r  =  tan'M—  ( Chirp  rate)) 
c 

where  X  and  c  ,  respectively,  denote  wavelength  and  the  velocity  of 
light.  Hence,  the  functional  form  for  the  projections  can  be  stated 
as  follows. 


P*U) 


_ 2r _ 

_i 

[Ar2  -  [AB  -C)  l2]  2 


A  »  sin2  if  +  cos*  t  ,  B  =  -il^2 +  cos2  if  ,  C  =  (l  -  -^)cosa  t  sin2  ♦ 
0  0  0a 

An  illustration  of  this  is  provided  in  figure  23b,  while  the 
tomographic  reconstruction  of  the  distribution,  for  a  given  set  of 
the  parmeters  r  and  0,  is  shown  in  figure  23c. 
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(a)  SYSTEM  GEOMETRY 


T«t« 


(b)  THE  RANGE ( R ) -RANGE  RATE ( R )  DISTRIBUTION  AND  A  SAMPLE 
PROJECTION  FOR  ANGLE  ^ 


(C)  RECONSTRUCTED  RANGE ( R ) -RANGE  RATE ( R )  DISTRIBUTION 


FIGURP  23  A  PRIMATIVE  EXAMPLE  OF  A  TOMOGRAPHIC 
RECONSTRUCTION  IN  RADAR,  IN  WHICH  A 
ROTATING  RING-LIKE  OBJECT  IS  THE 
TARGET  . 
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5.0  COMPARISON  OF  TWO  IMAGING  APPROACHES 


Gy  exploiting  coherent  integration  to  process  the  pulses  echoed 
back  from  some  target,  a  conventional  Doppler  processing  system 
(  figure  24  )  will  produce 
maximum  gain  in  signal  to  noise 
and,  theoretically,  can  achieve 
greater  Doppler  resolution  by 
adding  pulses  to  the  pulse  train 
that  is  transmitted.  In  spite  of 
these  seemingly  superior 
attributes,  the  investigation  of 
an  alternative  is  justified 
because  of  the  following 
situations . 


1.  Doppler  as  well  as  range  delay 
are  usually  not  constant. 

Normally,  there  is  steady 

migration  through  resolution  cells 
by  the  scattering  points, 

particulary  in  the  case  of 
rotating  targets.  As  a  result,  the 
effectiveness  of  the  coherent 
integration  diminishes  as  a  function  of  time. 

2 .  The  coherent  integration  breeds  multiple  pop-up  ambiguities  that 
are  characterized  in  the  literature 
as  a  "Bed  of  Nails"  (see  figure  25) 

Bounded  in  range  delay  by  the  period 
between  pulses  and  by  the  reciprocal  of 
the  period  with  respect  to  Doppler 
shift,  the  unambiguous  area  can  never 
be  enlarged.  Unless  the  area  to  be 
imaged  is  similarly  bounded,  the 
ambiguities  are  manifested  by  aliasing 
artifacts  superimposed  on  the  images. 


3 .  To  derive  the  maximum  performance 
that  is  offered  by  coherent 
integration,  a  rigid  schedule  of  pulse 
transmissions  is  required.  In  a 
multipurpose  system,  this  can  conflict 
with  other  transmissions  that  also  have 
to  be  scheduled. 

By  integrating  the  echoed  pulses  after 
detection,  tomographic -based  imaging 
offers  ambiguity- free  coverage  which  means  no  pop-up  ambiguities 
across  the  entire  domain  of  range  delay  and  Doppler  shift.  In 
addition,  the  pulse  transmissions  can  be  scheduled  independently 
which  lessens  the  possibility  of  transmission  conflicts  in  a 
multpurpose  system.  Also,  the  freedom  to  transmit  pulses  at  any 
time  can  be  utilized  to  achieve  better  snapshots  of  the  range- 
Doppler  distributions  by  speeding  up  the  transmission  of  pulses  as 
a  stop  action  measure. 

The  table  next  page  summarizes  made  above. 
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FIGURE  25  "BED  OF  NAILS" 
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COMPARISON  OF  TECHNIQUES 


Ducumoii 

TECHNIQUE 


STATUS 


TRANSMIT 


RECEIVE 


ISSUES! 

S/N. 

GAIN 


AMBIGUITY . 

CHARACTERISTIC 


EFFECTIVE. . 
INTEGRATION 
TIME 


PULSE . 

TRANSMISSIONS 


PHASE. . . . 
COHERENCE 


PULSE  DOPPLER  RADAR 

Conventional 

Wideband  Pulse  Train 

Pulse  Compression  & 
Coherent  Doppler 
Integration 

Maximum 

■BED  OF  NAILS* 

Limited  due  to 
target  dynamics 

Ridgid  pulse  schedule 

Intrapulse  stability 


CHIRP  DIVERSITY  RADAR 


New  Idea 


Diverse  Linear-FM 
Pulse  Suite 

Pulse  Compression  & 
Noncoherent  Tomographic 
Combining 


<  Maximum 


SUPERTHUMBTACK 


Snapshot  measurements 
offer  unlimited 
diversity  combining 


Independent  pulses 


Accurate  interpulse 
FM  rate 
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6 . 0  COMPUTATIONAL  FRAMEWORK 


To  reconstruct  a  two-dimensional  distribution  for  which 
measurements  have  provided  the  line-integral  projections  that 
constitute  its  Radon  Transform,  any  one  of  three  established 
tomographic  computational  procedures  can  be  employed.  These  consist 
of  the  (1)  the  Inverse  Radon  Transform,  (2)  computations  based  on 
the  Projection-Slice  Theorem,  and  (3)  an  algebraic-based 
reconstruction.  The  reconstructions  that  have  been  produced  during 
this  study  were  derived  by  filtering  the  projections,  then 
backprojecting  the  response.  Since  the  filtering  was  done  by 
appropriately  weighting  DFT  amplitudes,  these  tomographic 
reconstructions  were  actually  obtained  by  a  combination  of  the 
procedures  listed  above  as  (1)  and  (2)  .  This  computational 
framework  can  also  be  employed  for  systems  in  which  tomographic 
methods  are  combined  with  a  Chirp  diversity  radar  in  order  to  map 
radar  backseat ter  as  a  function  of  range  and  Doppler. 

Figure  26  presents  one  possible  system  configuration.  The  Chirp 
diversity  is  provided  in  this  system  during  a  sequence  of  linear  FM 
pulse  transmissions.  The  receiver  is  shown  with  a  common  pulse 
compression  network  which  is  preceded  by  a  mixing  circuit.  The  two 
inputs  consist  of  (1),  the  delayed  echo  with  linear  FM  rate  Vj,  and 
(2),  a  reference  pulse  with  a  small  increment,  Q,  in  the  Chirp 
rate.  The  latter  is  delayed  so  that  it  coincides  with  a  time  delay 
corresponding  to  a  reference  range. 

The  mixing  process  produces  an  output  pulse  which  is  approximately 
the  same  duration  as  the  input  but  with  the  linear  FM  rate  reduced 
to  £2.  The  purpose  is  twofold:  (1)  the  signal  bandwidth  is  reduced 
for  easier  processing,  and  (2)  the  FM  rate  is  normalized  to  enable 
the  pulse  compression  with  a  common  circuit. 
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Although  the  bandwidth  is  reduced  by  the  preceding  processing 
scheme,  resolution  in  range  is  not  sacrificed.  This  is  explained  by 
considering  the  output  time  delay  difference  for  two  different 
input  range  displacements.  Thus,  the  time  delay  difference  at  the 
output  is  scaled  by  the  same  factor  that  the  output  pulses  are 
widened  due  to  narrower  pulse  compression  bandwidths.  This  scale 
factor  equals  the  ratio  of  the  original  bandwidth  divided  by  the 
bandwidth  after  the  mixing.  By  maintaining  the  ratio  of  these 
bandwidths  with  every  pulse  transmission,  pulse  to  pulse  variations 
in  output  time  scale  are  avoided;  a  constant  time  scale  is  required 
for  inputs  supplied  to  the  tomographic  processor. 


The  preceding  "time  expansion*  system  yields  projections  over  range 
delay.  However,  analogous  projections  can  also  be  produced  over 
Doppler.  These  projection  options  are  illustrated  in  figure  27 
which  also  illustrates  the  linear-mapping  transformation  that 
arises  as  a  result  of  the  characteristic  coupling  of  range  delay 
and  Doppler  during  Chirp  pulse  compression  processing.  The 
coordinate  systems  involved  are:  (1)  range  (R)  and  velocity  (R) , 
and  (2),  range  delay  (X)  and  Doppler  offset  (0) .below. 


Expressed  in  matrix 
format,  the  linear 
relationship  between 
these  coordinates  is  as 
follows . 


T 
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Figure  27  displays  a  simple  scenario  consisting  of  three  isolated 
backscatter  sources  with  cross-sections  at,  o2,  and  o3.  This  is  is 
mapped,  simultaneously,  on  the  RR-plane  and  on  the  if-plane  which 
is  rotated  by: 


e  -  raw1  1 


-  {f.£} 


The  figure  can  be  interpreted  as  follows.  In  general,  the 
tomographic  projections  P9(x)  and  P9(0)  are  comprised  of  backscatter 
that  accumulates  in  windows  with  alignments  perpendicular  to  range 
delay  and  Doppler  offset,  respectively.  The  window  sizes  with 
respect  to  range  delay  and  Doppler,  respectively,  are 


(1) 


At 


1 

A  f 


(2) 


For  the  particular  case  that  is  illustrated,  both  projections 
consist  of  a  pair  of  dirac-like  responses,  while  for  the  specific 
windows  that  are  illustrated,  the  response  amplitudes  are 


(1)  Pe(  t)  =  ax  +  o2 

(2)  P«($)  •  ox  +  o3 


Some  other  Chirp  rate  and  corresponding  rotation  of  the  T^-plane 
will  yield  different  projections. 

As  described  above,  the  projections  on  x,  which  are  represented  by 
P©  ( X ) ,  are  available  as  outputs  from  a  system  such  as  in  figure  27. 
To  extract  the  analogous  projections  P9($),  it  is  necessary  in 
effect  to  corner-turn  the  pulse  compression  processing.  This  can  be 
done  by  (1)  ,  nulling  the  FM  increment  in  the  reference  Chirp  signal 
that  is  supplied  to  the  mixer  and  (2),  replacing  the  pulse 
compression  circuit  with  some  kind  of  spectrum  analyser  (  eg.,  FFT 
or  the  Chirp  Z  Transform)  .  The  schematic  for  the  revised  system  is 
provided  in  figure  28. 
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The  system  that  is  ultimately  chosen  to  image  backscatter 
distributed  over  range  and  Doppler  can  consist  of  either  one  of  the 
configurations  in  figures  26  and  28.  Another  option  is  a 
configuration  in  which  the  tomographic  computations  are  performed 
by  combining  the  projections  from  both  systems.  A  configuration 
such  as  this  is  attractive  since  it  offers  the  same  image  quality 
with  half  as  many  Chirp  transmissions. 
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7 . 0  CONCLUSION 


The  conceptual  system  for  conventional  mapping  of  radar  echo 
intensities  in  range  delay  and  Doppler  shift  coordinates  consists 
of  a  bank  of  matched  filters  whose  center  frequencies  are 
progrssively  offset  with  respect  to  the  carrier  frequency;  the 
equivalent  is  a  two-dimensional  array  of  correlators  for 
displacements  in  range  delay  as  well  as  Doppler  shift.  The  maps 
that  are  produced  are  never  exact,  however.  This  condition  is  a 
symptom  of  the  point  spreading  which  is  characterized  by  the  radar 
ambiguity  function.  The  point  spreading,  an  inescapable  law  of 
nature,  is  the  penalty  which  is  paid  because  the  view  of  objects  is 
always  through  windows  hence  the  information  available  for 
conventional  imaging  is  always  incomplete.  Lord  Rayleigh  recognized 
this  fact  and,  thereupon,  provided  his  well-known  resolution 
criterion  which  quantifies  our  inability  to  see  things  as  they 
really  are. 

Although  optimization  is  theoretically  available  to  the  designer 
who  can  specify  idealized  waveforms,  he  ultimately  must  accommodate 
certain  constraints  as  well  as  conflicting  system  requirements  that 
limit  how  much  can  be  achieved  in  practice. 

The  practical  solution,  up  until  now,  has  been  to  combine  the 
echoes  returned  to  radar  corresponding  to  a  sequence  of  coherent 
pulse  transmissions.  These  transmissions  are  often  encoded  by 
internal  modulation  so  as  to  enable  pulse  compression  when  the 
echoes  are  received.  The  system  is  referred  as  pulse  Doppler,  and 
the  corresponding  point  spread  function  is  commonly  characterized 
as  a  "bed  of  nails". 

This  pulse  Doppler  solution  normally  only  provides  access  to  the 
central  region  of  the  ambiguity  plane.  The  boundaries  are 
established  by  the  pulse  repetition  frequency  ( PRF)  in  one 
direction  and  by  the  period  between  the  pulses  in  the  other 
dimension.  Any  echoe  intensity  from  a  source  outside  this  region  is 
folded  back  onto  the  central  region.  This  fold-back  intensty  cannot 
be  distinguished  from  true  distribution  of  intensity  and  is 
commonly  termed  "aliasing". 

A  serious  attempt  is  usually  made  during  system  design  to  limit  the 
amount  of  aliasing  that  can  occur.  Idealistically,  the  central 
region  is  made  sufficiently  large  so  that  no  intensity  distribution 
extends  beyond  the  boundaries.  This  goal  is  often  elusive,  however. 
The  reason  being  the  reciprocal  coupling  that  exists  between  the 
PRF  and  the  pulse  period,  which  bars  simultaneous  enlargement  of 
both  dimensions  of  the  central  region.  For  many  situations, 
therefore,  the  mapping  accuracy  will  depend  on  the  final  compromise 
that  is  acceptable  --  all  things  considered. 

The  likelihood  for  such  situations  as  above  arising  during  the 
design  of  the  conventional  system  can  be  expected  to  grow  in  the 
future.  This  will  be  driven  by  a  requirement  for  higher  carrier 
frequencies  to  accommodate  the  wider  pulse  compression  bandwidths 
that  will  be  necessary  to  meet  expected  demands  for  greater  range 
resolution.  For  a  range  resolution  of  less  than  one  foot,  it  will 
be  necessary  to  employ  millimeter  wave  transmissions  to  mitigate 
the  spectrum  distortion  at  the  low  end  due  to  fold-over.  The 
unambiguous  detection  area  will  shrink,  as  a  result,  due  to  greater 
Doppler-velocity  'sensitivity .  A  compounding  of  this  dilemma  can 
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occur  if  a  greater  unambiguous  detection  range  is  also  demanded 
which  will  require  a  longer  pulse  period.  The  probability  that 
there  will  be  faster  targets  in  the  future,  further  intensifies  the 
likelihood  their  images  will  be  ambiguous. 

By  integrating  the  echoed  pulses  after  detection,  a  Chirp  diversity 
radar  incorporating  tomographic  methods  posseses  none  of  •these  -77^2 
symptoms  associated  with  conventional  pulse  Doppler  techniques.  The 
point  spreading  is  a  thumbtack- like  function. 

Thus,  range-Doppler  coverage  is  unlimited  since  there  are  no 
ambiguity  pop-ups  to  restrict  this  coverageg.  In  addition,  it 
appears  that  the  mainlobe  width  depends  exclusively  on  bandwidth. 

In  contrast,  the  classic  thumbtack  ambiguity  function  shown — i*u 
— 2_  is  affected  both  by  bandwidth  and,  inversely,  by  the 
integration  time.  More  important  than  the  preceding  discovery  is 
the  property  that  ambiguity  function  is  not  constrained  by  radar 
uncertainty  principles.  Hence,  "the  sands  of  the  beach  *  is  not  a 
suitable  analogy  to  descibe  this  smbiguity  function.  It  is 
concluded  that  this  superthumbtack  characteristic  offers  superior 
resolution  for  an  accurate  interpretation  of  radar  backscatter  with 
respect  to  detection,  parameter  estimates,  and  imaging. 

An  added  benefit  by  integrating  the  echoed  pulses  after  detection 
is  the  freedom  to  schedule  pulse  transmissions  independently.  The 
advantages  gained  as  a  result,  over  conventional  imaging 
techniques,  are  the  following.  First,  the  freedom  to  tansmit  pulses 
at  any  time  enables  accelerated  measurements  which  can  yield 
images  with  stop-action  sharpness.  Second,  since  the  pulse 
transmissions  are  not  constrained  to  a  ridgid  PRF  schedule,  they 
can  be  interleaved  with  other  transmissions  in  a  multifunction 
radar  so  as  to  utilize  available  energy  in  the  most  efficient 
manner . 


The  radar  applications  in  which  the  Chirp  diversity  combined  with 
tomography  can  be  employed  are  many.  In  general,  this  covers  all  of 
the  categories  where  range  and  Doppler  information  is  utilized. 
Some  examples  are  ICBM  wake  classification,  surveillance  and 
tracking,  weather  mapping  and  microburst  detection,  traffic 
control,  ionospheric  sounding,  radar  astronomy,  and  imaging 
rotating  targets.  The  value  of  this  technique  will  become 
significant  when  it  is  no  longer  possible  to  exploit  coherent 
techniques  effectively.  The  environment  for  this  to  unfold  will 
occur  perhaps  in  a  future  SDI  scenario  or  possibly  in  probing  outer 
space  for  menacing  objects.  Laser  radars  operating  at  optical 
freqencies  are  particulary  vulnerable  regarding  ambiguous  Doppler, 
while  potentially  high  target  speeds  compared  acoustic  propogation 
velocity  makes  sonar  simarly  vulnerable.  It  is  predicted  that  this 
application  of  tomography  will  eventually  become  the  principal 


radar  approach  for 
and  Doppler  shift. 


lapping  radar  backscatter  as  a  function  range 
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APPENDIX  A 


A  TEST  AND  DEMONSTRATION  CONFIGURATION 


I.  RAYTHEON  FACILITIES 


Raytheon  Company  at  Sudbury  and  Wayland,  Massachusetts  currently 
maintains  several  facilities  for  testing  microwave,  millimeter- 
wave  ,  and  laser  radars.  They  are  conveniently  placed  in 
laboratories  inside  the  plant  as  well  as  on  surrounding  grounds, 
and  are  regularly  utilized  for  research  and  development  in 
connection  with  new  antenna  system  designs  design  as  well  as  for 
predicting  the  radar  cross-section  of  different  targets. 

Effective  utilization  of  these  facilities  is  ensured  by  resident 
teams  consisting  of  qualified  mirowave  engineers  and  technicians. 
In  addition,  the  instrumentation  equipment  permanently  installed 
at  each  facility,  consisting  of  sophisticated  microwave  measuring 
devices  as  well  as  computers,  also  ensures  precise  measurements 
and  a  rapid  reduction  the  data  to  useful  information. 

Any  of  these  facilities  offers  a  testing  environment  that  is 
suitable  for  controlled  experiments  which  can  demonstrate  the 
capability  tomographic  imaging.  But  there  may  be  similar 
facilities  not  owned  by  Raytheon  Company  which  will  require  less 
costly  developments  in  order  to  implement  a  test  and 
demonstration  program.  An  investigation  of  the  possibility  is 
certainly  worth  considering. 

II.  SIGNAL  PROCESSING  FUNCTIONS 


The  functional  diagram  for  a  possible  test  configuration  is 
illustrated  in  figure  Aa.  Functions  are  divided  between  tasks 
performed  with  a  computer  and  various  signal  processing  tasks. 
The  former  consist  of 
tomographic 
computations  and  also 
computations  which  will 
produce  a  suite  of 
linear  FM  waveforms 
with  different  chirp 
rates.  By  incorporating 
a  provision  for 
transporting  digitized 
signals  via  magnetic 
tape,  these 
computations  can  be 
performed  off -site. 

This  will  enable 
ongoing  analyses  to 
continue  at  the 
workstation,  without  an 
int erupt ion  in  order  to 
relocate  the  computer 
to  the  test  site. 

To  buffer  the  digitized 
signals,  a  RAM  memory 
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will  be  installed  at  the  test  site  with  sufficient  capacity  to 
store  input -output  data  for  one  radar  transmission.  It  will  be 
necessary  for  this  digital  memory  to  provide  a  buffer  for  high 
speed  (  500  MHZ  bandwidth  )  output  data  as  well  as  for  relatively 
narrow  bandwidth  data  of  less  than  500  KHZ. 

The  wideband  output  data  will  consist  of  (12  bit)  I  and  Q  digital 
samples  that  are  transported  from  the  computer  via  magnetic  tape. 
Representing  one  chirp  modulation  pulse,  this  data  will  be  output 
to  a  D/A,  then  channeled  to  a  single  sideband  modulator  and, 
subsequently,  supplied  as  as  input  to  a  (60  MHZ  )  microwave 
transmitter. 

The  narrow  bandwidth  input  data  will  consist  digitized  (12  bit)  I 
and  Q  samples  representing  'stretch'  processed  radar  echoes  for  a 
single  chirp  transmission.  (  The  'stretch  processing  is  a 
familiar  technique  whereby  the  echoes  are  dechirped  by  mixing 
them  with  a  an  undelayed  replica  of  the  signal  transmitted  ) . 

III.  TRANSMIT/RECEIVE  CONFIGURATION 


Separate  transmit  and  receive  horns  are  required.  This 
requirement  arises  because  their  signals  overlap  and  must 
therefore  be  isolated.  The  horns  will  be  elevated  to  an  altitude 
of  60  feet  and  made  to  face  downward  to  a  rotating  platform  on 
the  ground  at  approximately  134  feet  radar  range  (  figure  A2  )  .  A 
similar  installation  has  been  employed  at  the  Environmental 
Research  Institute  of  Michigan  for  studies  concerning  SAR  and 
ISAR  signal  processing  [*] .  The  specification  here  is  that  the 
platform  have  20  feet  diameter  and  be  able  to  rotate  at  a  maximum 
rate  of  30  RPM.  This  will  provide  a  maximum  Doppler  shift  of 
several  kilohertz  after  compensation  for  the  oblique  intercept 
angle.  (  The  maximum  Doppler  shift  should  be  sufficient  for  the 
experiments  that  are  to  be  performed.)  Adequate  platform 
stabilization  is  required  with  respect  to  speed  as  well  as 
transverse  motion,  which  will  limit  spurious  Doppler  spreading  to 
a  maximum  of  100  HZ. 

IV.  HORN  DESIGN 


The  design  of  the  horns  for  transmitting  and  receiving  will 
provide  radiation  beams  that  will  cover  the  entire  surface  area 
of  the  rotating  platform.  Approximate  horn  dimensions  required 
for  this  at  0.5  cm  (60  GHZ)  wavelength  are  as  follows: 

Length  =15  cm 

H  Plane  Aperture  =5  cm 

E  Plane  Aperture  =  4.34  cm 

V.  TRANSMITTER  POWER 


To  provide  20  DB  echoes  relative  to  noise,  for  a  6  square 
centimeter  RCS  mounted  on  the  rotating  platform  target  requires 
15  mw  peak  power.  This  is  based  on  27  DB  gain  computed  for  the 
horns  and  assumes  41.7  DB  pulse  compression  gain  is  obtained  by 
transmitting  a  30  us  pulse  that  is  frequency  swept  over  a 
bandwidth  of  500  MHZ.  This  is  a  maximum  requirement,  since  the 
peak  power  requirement  will  drop,  assuming  constant  bandwidth,  as 
the  pulse  duration  is  increased  to  acquire  lower  FM  slopes  for 
full  angle  coverage  regarding  tomographic  projection  aspects. 
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FIGURE  Aa 


a.  TARGET  MODEL 


b.  WAVELENGTH  _ 

c.  XMT  POWER  . 

d.  MODULATION  _ 

e.  SWEPT  BANDWIDTH 

f.  FM  RATE  . 


f.  SIGNAL  PROCESSING  . 

h.  TOMOGRAPHIC  ALGORITHM . 

i.  PLATFORM:  (1)  SPEED  . 

(2)  STABILITY.. 


(1)  Point  Target  (6  cm*  RCS) 

(2)  Multiple  Point  Targets 
(various  RCS  per  target) 

(3)  Ring  Target 

Millimeter-wave  (60-100  GHZ) 
15mw  for  20  DB  S/N 
Chirp  Diversity 
500  MHZ  (nominal) 

.+/-  megahertz  per  ps 
Long  CW  Pulse 
0.10) 


6.50 » 

Short  CW  Pulse 

"Stretch"  Pulse  Compression 
Filtered  Backproject ion 
1/2  rps 
100  HZ 


[*]  J.L.  WALKER, Range -Doppler  Imaging  of  Rotating  Objects, 
Dissertation  for  D.  Phil.,  U.  of  Michigan,  (1974) 
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